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Spin effects in Kapitza-Dirac scattering at light with elliptical polarization 


Rico Erhard and Heiko Bauk^ 

Max-Planck-lnstitut fur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 

The Kapitza-Dirac effect, which refers to electron scattering at standing light waves, is studied in the Bragg 
regime with counterpropagating elliptically polarized electromagnetic waves with the same intensity, wavelength, 
and degree of polarization for two different setups. In the first setup, where the electric field components of the 
counterpropagating waves have the same sense of rotation, we find distinct spin effects. The spins of the scattered 
electrons and of the nonscattered electrons, respectively, precess with a frequency that is of the order of the 
Bragg-reflection Rabi frequency. When the electric-field components of the counterpropagating waves have an 
opposite sense of rotation, which is the second considered setup, the standing wave has linear polarization, and no 
spin effects can be observed. Our results are based on numerical solutions of the time-dependent Dirac equation 
and the analytical solution of a relativistic Pauli equation, which accounts for the leading relativistic effects. 

PACS numbers: 03.65.Pm, 31.15.aj 


1. Introduction 

Kapitza-Dirac scattering lEEl, which is the scattering of mas¬ 
sive particles on a standing light wave, emerges as a remarkable 
consequence of the wave-particle duality of quantum objects. 
When Kapitza and Dirac first predicted electron scattering at a 
standing light wave in 1933, an experimental realization was 
out of reach. The advent of high-intensity laser facilities, how¬ 
ever, rendered such an experimental realization possible. The 
first clear observations of the Kapitza-Dirac effect in the so- 
called Bragg regime were achieved by employing atoms El SI 
in 1986 Q and in the so-called diffraction regime |0| in 1988. 
Recently, strong-field Kapitza-Dirac scattering of neutral atoms 
was realized through frustrated tunneling ionization leading to 
an extreme acceleration of the employed helium atoms Gl. The 
Kapitza-Dirac effect with electrons was first achieved in 1988 
(8| in the diffraction regime. The scattering of electrons in the 
Bragg regime was demonstrated experimentally in 2001 Elllol, 
which comes closest to the diffraction process as proposed orig¬ 
inally by Kapitza and Dirac m. These experiments stimulated 
also renewed theoretical interest in this effect unHni and gen¬ 
eralizations thereof For example, in Refs. nmia two-particle 
Kapitza-Dirac scattering was studied, and also Kapitza-Dirac 
scattering of electrons from a bichromatic standing light wave 
with linear polarization IfT^ or circular polarization IfTTlI has 
been considered recently. In dielectric media, electrons may 
also diffract from traveling light waves im. 

As the electron couples to electromagnetic fields not only 
via its charge but also via its spin degree of freedom, it is 
natural to ask if electron-spin dynamics may be observed in 
the Kapitza-Dirac effect ||9l. Freimund and Batelaan found 
a vanishingly small spin-flip probability in the investigated 
laser setup and parameter regime by simulating the Bargmann- 
Michel-Telegdi equations ini. Another derivation ll^ solved 
the Pauli equation perturbatively via second-quantization meth¬ 
ods in the diffraction regime but only found tiny spin effects 
for a setup with linearly polarized light waves. Spin effects. 
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however, can be observed in the Kapitza-Dirac effect if the 
electron enters the laser field with a certain relativistic initial 
momentum, which must fulfill a spin-flip resonance condition 
l|2l]|22l, or in specific laser setups with circular polarization 

na. 


As a result of emerging novel high-intensity light sources 
such as the Extreme Light Infrastructure Ultra High Field Fa¬ 
cility, which envisage to provide field intensities in excess of 
I0^®W/cm^ and field frequencies in the x-ray domain ll23l . 
spin effects in relativistic light-matter interaction became of 
wide interest ll24lj^ . In particular, the interaction between the 
electron’s spin and the degree of elliptical polarization of an 
external strong electromagnetic field has been considered. For 
example, it was recently shown that electrons scattering non- 
resonantly at standing light waves of counterpropagating plane 
waves with elliptical polarization can exhibit spin precession 
Il29l|30l. Strong spin effects have also been found in electron- 
positron pair production in strong electromagnetic fields with 
elliptical polarization RTl [3^ as well as in strong-field ion¬ 
ization via electromagnetic fields with elliptical polarization 

lESllMl. 


In this paper we will study the Kapitza-Dirac effect with 
electrons in standing light waves of counterpropagating plane 
waves with elliptical polarization and especially study pos¬ 
sible spin effects. This article is organized as follows; In 
Sec. we will specify the considered configurations of the 
time-dependent electromagnetic fields. The effect of the fast 
oscillating fields can be described by time-independent pon- 
deromotive potentials, which will be derived in Sec.[^ The 
resulting time evolution by the ponderomotive potentials will be 
determined analytically and compared to a numerical solution 
of the time-dependent Dirac equation in Sec. A summary 
and our conclusions follow in Sec.|5] 
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2. Electrons in counterpropagating 
electromagnetic waves with 
elliptical polarization 

Electron scattering at a standing light wave, which is composed 
of two counterpropagating laser waves with equal wavelength 
A, equal electric field amplitude E, and the same degree of 
ellipticity, will be considered. Depending on the relative sense 
of rotation of the electric-field component of the two plane-wave 
fields traveling along the x axis, the electric- and magnetic-field 
components are given by 
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respectively, where c denotes the speed of light, k = InjA is 
the wave number, and to = IncjA. The parameter rj e (- 7 r, 7 r] 
determines the degree of ellipticity, with rj = 0 and rj = ti cor¬ 
responding to linear polarization and 77 = ±7r/2 corresponding 
to circular polarization. The two electromagnetic waves given 
by ([T^ and ( [Tbl l rotate in parallel (corotating setup), whereas 
the two electromagnetic waves given by ( | 2 al l and ( | 2 b| ) rotate 
opposite to each other (antirotating setup); see Fig.[ 2 for an illus¬ 
tration. The total electromagnetic fields and the corresponding 
vector potentials of these setups follow as 
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FIG. 1: (Color online) Schematic illustration of the two considered 
field configurations. Two electromagnetic plane waves with elliptical 
polarization traveling into opposite directions are superimposed. The 
sense of rotation of the electric field vector may be (top) equal (coro¬ 
tating setup) or (bottom) opposite (antirotating setup). The sense of 
rotation and the propagation direction are indicated by black arrows, 
red solid arrows represent electric-field components, and blue dashed 
arrows represent magnetic-field components. 


A^^(x, t ) 


2E 

to 


sin(mf) 


0 

cos{kx) 
^cos(^x + 77)^ 


(4c) 


for the antirotating setup, respectively. 

The quanrnm-mechanical evolution of an electron of mass 
777 and charge q = -e in the laser field with the vector potential 
( [^ or ( |4c] l is governed by the quasi-one-dimensional Dirac 
equation 


i/ 7 'i'(x, t) = ^-cihox- - cqw{t)a ■ A(x, t) + mc^/3^ 'I'(x, t) 


( 5 ) 


with the Dirac matrices a = {ax,ay,a^)^ and (5 llTSlI^ and the 
reduced Planck constant h. This form of the Dirac equation can 
be justified by starting from the fully three-dimensional Dirac 
equation and noting that the canonical momentum is preserved 
in the y and in z directions due to the special form of the vector 
potentials ( |^ and ( |4cl i. Thus, the three-dimensional Dirac 
equation separates into three one-dimensional equations, of 
which one is given by and the remaining two describe free 
motion in the y and in z directions. In (|^ the window function 
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has been introduced, which allows for a smooth turn-on and 
turn-off of the laser field. 

In Sec. 1^ we will present numerical solutions of the time- 
dependent Dirac equation 0 which are obtained solving this 
equation on a one-dimensional grid of size A with periodic 
boundary conditions by a Fourier split-operator method 037II38I . 




























3 


Due to the periodic boundary condition the wave function can 
be expanded into a discrete set of momentum eigenstates 

(7) 

which are simultaneous eigenstates of the free Dirac Hamil¬ 
tonian, the X component of the canonical momentum op¬ 
erator (with eigenvalue nkfi), and the z component of the 
Foldy-Wouthuysen spin operator |[39l. The upper index y e 
{+t, -T, +i, -i} in ij/Z (jc) indicates the spin state, and the sign 
of the energy eigenvalue and uZ is defined as 



with;^^i = (1,0)^ and^^ = (0,1)^ and the relativistic energy- 
momentum relation 

£n = \/{nic^Y + {nckhy . (9) 

In the so-called Bragg regime of the Kapitza-Dirac effect cmii 
the electron has to meet the resonance condition, which requires 
that the electron’s momentum in the propagation direction of 
the light field equals ±hk. Thus, we will employ as a 

canonical initial quantum state. 


3. Ponderomotive potentials 


with the magnetic vector potential A{r,t) and B{r,t) = 
V X A(r, f), E{r, t) = -A{r, t), and the vector of Pauli matri¬ 
ces <T = {cr^,cry,cr^y. Taking into account that the electron’s 
kinematic momentum remains small compared to me and that 
/I » /Ic = hl{mc) allows us to neglect various terms in Eq. ( [T0| ), 
and the specific form of the vector potentials and ( |4c| l leads 
finally to 


m(^j) . (- - *<T ■ B(x.,) 

^ ^ \ 2m dx^ 2m ^ ^ 2m ^ ' 


qh 

4m^c^ 


cr ■ {E{x, t) X (-i/iVx - qA{x, ())) 1 'I'(x, f), (11) 


where Vx denotes the differential operator Vx = (<9x,0,0)^. 
This is the nonrelativistic Pauli equation, which is amended 
by a relativistic correction due to spin-orbit coupling 
(~ £■ X (-i/zVx)) and the coupling of the electron’s spin to the 
photonic spin density of the electromagnetic wave (- £ x A) 
Il29l[30l . It incorporates the dominating terms of ( [TOl i that ac¬ 
count for the electron’s spatial motion as well as for the spin 
dynamics. 

As shown in ina for the case of linear polarization, the 
Kapitza-Dirac effect can be modeled by a time-independent 
Hamiltonian with suitable ponderomotive potentials. Thus, it 
appears attractive to replace the time-dependent Hamiltonian 
in Eq. ( [TT] i by a time-independent Hamiltonian, which we will 
derive in the following by means of a Magnus expansion HTl 
\42\ . The Magnus expansion gives an exponential representation 
of the solution of a first-order homogeneous linear equation of 
the type 


The Bragg regime of the Kapitza-Dirac effect is realized for 
sufficiently weak electromagnetic fields, more precisely, if the 
energy that is gained in a single laser cycle, which is of the 
order c\q\E/ai, is small compared to the energy of a single 
photon, which is ha>. Furthermore, we will assume that the 
electron’s kinetic energy remains small compared to its rest 
mass energy mc^ for all times. Thus, c\q\E/a> « hat < mc^. As 
a consequence of the Bragg condition, nonrelativistic initial mo¬ 
mentum implies T » Tc = (me), where Ac is the electron’s 

reduced Compton wavelength. The nonrelativistic electron 
momentum justifies the application of a weakly relativistic the¬ 
ory. The weakly relativistic limit of the Dirac equation may 
be reached via a Foldy-Wouthuysen transformation Il39ll40l . 
which yields, for a now two-component wave function 'I'(r, f), 


f(t)=H{t)f(t) (12) 

in the form 

/(f) = exp(f7i(f) + U 2 {t) + ... )/(0). (13) 

The main advantage of the Magnus expansion is that each 
truncated series ( |T3] l of any order preserves the unitary character 
of the quantum-mechanical time evolution. Explicitly, the first 
two terms of the series o are 

f7i(f)= r H{h)dt, (14) 
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The term U 2 {t) and higher-order corrections account for the 
fact that H{t) may not commute with itself at different times f 
and in this way implement time ordering. 

Applying the Magnus expansion in second order to ( [TT] ) 
with the vector potential of the corotating case ( [^ and the 
corresponding electromagnetic fields yields 
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where we have neglected in the exponent terms which do not grow linearly with time and/or are small (compared to the included 
terms) in the Bragg regime, i. e., \q\El{a)mc) « 1 with A» Ac = Til{mc). Equation ( [T6] l represents the time-evolution operator of 
the equation 


i/!'i'(x, f) = ^ , cos^{kx) - ^ sin 77 (sin^(/:x) - cos^{kx)) 'I'(x, f), (17) 

\ 2m dx^ map trPPu ) 

which is the searched-for evolution equation with a time-independent Hamiltonian. It equals the Hamiltonian for the case of linear 
polarization lfT3l plus a term which accounts for the spin dynamics due to the ellipticity of the electromagnetic field. Thus, we 
found polarization-dependent ponderomotive forces similar to the situation in Refs. MM- Note that the spin term originates in 
equal magnitude from both the Zeeman term --- B -cr and the relativistic correction ~ (E x A) ■ cr (see also |l29]|30l). 

Similarly, we can apply the Magnus expansion to o with the vector potential of the antirotating case ( |4cl i and the corresponding 
electromagnetic fields, which yields now 
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where we have again neglected in the exponent terms which do not grow linearly with time and/or are small in the Bragg regime. 
Equation ([TF| represents the time-evolution operator of the equation 
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This can be further simplified by a gauge transform, which 
removes the additive constant, and a suitable shift of the coor¬ 
dinate system to 
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In contrast to ( [T7] i, the Hamiltonian in ( [20l i involves no spin¬ 
coupling terms, which is a consequence of £ x A = 0 and 
the fact that the Zeeman interaction is zero on average (as it 
commutes with itself at different times) in the antirotating case. 


4. Analytical solution of the 
relativistic Pauli equation and 
numerical results 

In the following, we will derive approximate analytical solu¬ 
tions of the relativistic Pauli equations ( [T7| l and ( |20| l. Eor this 
purpose the wave function 'I'(x, t) is expanded into plane waves 
as 

( 21 ) 

with the time-dependent coefficients cl(t) and c;J(f), which 
will be conveniently combined into the pair 

Cr,(t) = (cl(t)’4(t)y ■ (22) 


The basis functions that are employed in ( [ST) represent states 
with a well-defined momentum and spin orientation. They 
form a system of orthonormal functions on x € [0,4]. Thus, 
the evolution equation for c^(f), the relativistic Pauli equa¬ 
tion in momentum space, can be derived by plugging the 
ansatz ( | 2 T] ) into ( [TTI i or ( | 20 l i and taking the scalar product with 
^Jkl{2n)x^ e‘" We will ignore possible turn-on and turn¬ 
off phases in the following, i. e., AT = 0, if not otherwise 
identified. 


4.1. Corotating fields 


Using the expansion ( |2T| l, the relativistic Pauli equation for 
corotating fields (fTT]) is given in momentum space by 
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Odd and even modes are decoupled in Eq. ( |23] l. Eurthermore, 
the relation 
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holds for the coefficients in Eq. ( |2^ because of the Bragg- 
regime condition c\q\E/a) « hoj < mp-. Thus, it is justified to 
truncate the equation system, which yields for the odd modes 
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by introducing the frequencies 
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and = Q 3 sin? 7 . Note that an additive constant 202 has 
not been included on the diagonal of the coefficient matrix in 
Eq. ( |25| l, which can be justihed by a suitable gauge transform. 
The eigenvalues and eigenvectors of this coefficient matrix can 
be easily calculated by a computer algebra system giving some 
rather intricate expressions. In the parameter regime, where 
( |24l i is fulhlled, the eigenvectors of the coefficient matrix can 
be approximated by the column vectors of the constant matrix 




V2 »3 »4 »5 »6 


Vi Vs) = 


1 

2 


/O 

0 

1 

1 

1 

1 

0 


\0 


0 0 0 

0 0 0 

1 1 1 

-1 -1 1 

-1 1 -1 

1 -1 -1 

0 0 0 

0 0 0 


1 1 

1 1 

0 0 

0 0 

0 0 

0 0 

1 -1 

1 1 


1 

1 

0 

0 

0 

0 

1 

-1 


1 

1 

0 

0 

0 

0 

-1 

- 1 / 


(27) 


and the corresponding exact eigenvalues are 
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Choosing c4(0) = 1 and c^(0) = 0 otherwise as the initial 
condition and approximating the eigenvectors of the coefficient 
matrix in Eq. @by ( |27l l yields the time-dependent solution 
to ( 1 ^, 


43 (f) =0, (29a) 
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t/l 

Note that although the amplitudes (f) are zero, also in¬ 
cluding modes n = ±3 in the truncated system ( |25| ) is crucial 
to describing the Kapitza-Dirac effect properly for corotating 

t/l 

helds. The amplitudes Cj /3 (f) vanish only because the approxi¬ 
mation ( |27l i has been applied. Employing the exact eigenvectors 
instead would yield small nonzero amplitudes. 

The evolution of the probability |4i(f)P is with (|29]l 
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Considering the regime ( |24l l, the following approximations 
hold: 
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Employing these approximations in Eq. ( [30l l, some algebraic 
transformations yield 
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A similar calculation can be carried out for |c4(f)p, |c}(f)4 
and |c|(f)p, which hnally results in 
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FIG. 2: (Color online) Time evolution of the probabilities to find the electron, which interacts with a standing light wave of the corotating setup, in 
a particular quantum state. Green (gray) lines show, from top to bottom, the analytical results p2a[ ( to p2d) based on the relativistic Pauli equation 
CD- Black lines represent the corresponding numerical results based on the Dirac equation (|^, i. e., the probabilities to find the electron in the 
states \ and respectively, after the electromagnetic field has been turned oif. Parameters are E = 400 a.u. = 2.06 x 10'“* V/m 

(corresponding to an intensity of 1.12 x 10^^ W/cm^), 4 = 3 a.u. = 0.159nm, AT = lOnjcj, and rj = njl, corresponding to circular polarization. 
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Thus, the occupation probabilities |cl[(f)p, |0j(f)p, |c|(f)p, 
and |c|(f)p show oscillatory behavior on different time scales. 
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The fast oscillation’s amplitude is modulated by the inter- 
mediate-frequency oscillation, and the fast oscillation’s am¬ 
plitude and the intermediate oscillation’s amplitude are modu¬ 
lated by the slow-frequency oscillation. The superposition of 
these three amplitude-modulated oscillations creates a nontriv¬ 
ial temporal behavior of the occupation probabilities (see also 
Fig.|2j. 

For ?7 = 0, i. e., £^3 = 0, Eq. ( [3^ simplifies to the known result 
for the Kapitza-Dirac effect with linearly polarized light ina. 


For nonzero ellipticity, however, quantum transitions occur not 
only between states with different momenta but also between 
states with different spin orientations. As a consequence, the 
standing wave’s ellipticity induces a beating behavior of the 
non-spin-resolved probabilities, i. e., 

kli(f)P + kii( 0 P ~ cos( 2 Q 20 cos( 2 £ 2 (f)) 

= i cos2((Q2 + + 5 cos^((Q 2 - , (36) 


kl(f)p + ki(0p ~ “ cos( 2 El 2 t) cos( 2 £ 2 (f)) 

= i sin^((02 + + 5 sin^((Q 2 - ■ (37) 

With ( |29) the spin expectation value 
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FIG. 3: (Color online) Time evolution of the spin of the two diffraction modes. Green (gray) lines show, from top to bottom, the analytical 
results ( |42[> and j43) based on the relativistic Pauli equation jl7) . Black lines show the corresponding numerical results based on the Dirac 
equation Parameters are as in Fig.|^ 
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FIG. 4: (Color online) Time evolution of the electron’s spin. Green 
(gray) line shows the analytical result ( |40) based on the relativistic 
Pauli equation ( |17) . Black line shows the corresponding numerical 
data based on the Dirac equation Parameters are as in Fig.|^ 


Thus, in contrast to the Kapitza-Dirac effect with linearly po¬ 
larized light, the electron’s spin precesses in a standing wave 
of elliptical polarization. If the spin is conditioned to a specific 
mode, let us say n = - 1 , one finds 




fi kli(0P-|cii(0P 

Ti (cos(2f22f) + cos(2f23f)) 

2 co&{2il2t) co&{2D!-^t) + 1 


(42) 


and, similarly, for n = 1 




it kKorkUop 

2kU0k + ki(0k 

h (cos(2Q20 - cos{ 2 Q.'^t)) 

2 cos(2f220“ 1 


(43) 


The analytical results, which have been obtained by deriv¬ 
ing an approximate solution to the relativistic Pauli equation 
shall be compared to numerical solutions of the time- 
dependent Dirac equation (|^ with space- and time-dependent 
electromagnetic potentials. Figure presents the probabili¬ 
ties to find the electron, which interacts with a standing light 
wave of the corotating setup, in one of the quantum states 


and and compares this to the analytical predic¬ 
tions ( |32a| ) to ( |32d| ). Even for rather long interaction times 
the analytical results ( |32a| ) to ( |32d| l match the fully relativistic 
dynamics as predicted by the Dirac equation. For the time evo¬ 
lution of the expectation value of the spin of the two diffraction 
modes, (s_i(f)) and (sj (f)), as well as for the time evolution of 
the expectation value of the spin (s(f)) we also find excellent 
agreement between the Dirac theory and ( |42l l, ( |43| l, and ( |40l l, 
respectively (see Figs. an dQ. The spin expectation values 
@, and ( |4^ oscillate on different time scales. Because 
Q 2 > shortest time scale of the total spin’s oscillation, 

which is determined by the frequency 2 Dk longer than the 
time scale of the oscillation of ( |42] i and ( |43] l, which is deter¬ 
mined by 2 Q 2 - Note that the electron’s spin precesses faster in 
the Kapitza-Dirac effect compared to nonresonant scattering at 
the same field configuration, where the spin precesses with the 
frequency D.2^3l{2-i^i) (see Refs. Il29l[^ and the ). 

Although we employed in the presented numerical examples 
for the Kapitza-Dirac dynamics ultra strong laser parameters 
in the soft x-ray regime, the observed electron-spin dynamics 
is not a strong-held effect. Spin hips in the diffracted beam 
may be observed also for less intense laser helds, but the spin- 
precession frequency becomes small in this regime. In this case 
the electron must be trapped for many laser cycles in the laser 
held’s focus to observe a full spin hip. The detectability of spin 
effects depends on the sensitivity of the spin measurement and 
for how long the electron interacts with the electromagnetic 
held. If one requires that the electron shall have reached full 
spin hip after Af laser periods and the laser parameters remain 
in the Bragg regime, then the required electric held strength is 
bounded from below and from above as 




< E « 


\q\c 


(44) 


Flere the lower bound follows from Af > a>I (402), with 2 O 2 
determining the shorter time scale of the spin dynamics of 
the diffracted electron beam. Note that the condition ( |44| ) is 
much less restrictive than the corresponding relation for the 
nonresonant scattering as presented in Ref. |30l. 

According to Eq. ( |40l ) the electron’s spin oscillates with 
the frequency 2 O 3 sin ij, which is modulated by an oscillation 
with the frequency O 2 O 3 sin?;/(20i). This allows us to test 
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FIG. 5: (Color online) Spin-oscillation frequencies for different de¬ 
grees of ellipticity 7 / as obtained by a numerical fit procedure to data of 
numerical solutions of the Dirac equation compared to the predictions 
by the relativistic Pauli equation ([TT). Parameters are as in Fig.|g 


the effect of the light’s ellipticity 77 on the Kapitza-Dirac ef¬ 
fect. For this purpose the numerical data for the spin expecta¬ 
tion value as a function of the interaction time were htted to 
/icos(Wflf) cos(aJif)/2 with the fit parameters a>a and ajb- In 
Fig.j^the numerical values of these ht parameters are compared 
to 2Q.J, and 02 ^ 3 /(20i). The numerical data show a clear sin rj 
dependency, as predicted by our analytical considerations. 


4.2. Antirotating fields 


For the setup with antirotating helds, which is described by 
( | 20 l i, the ansatz (| 2 T| gives the following set of equations in 
momentum space: 


ihcn{t) 


2m 


c„{t) + 


2 1^2 

q E cos Tj 
2k^mc^ 


{c„-2{t)+2Cn{t)+C n+2 it))- 


(45) 

Similar to the corotating case, only odd and even modes couple 
to each other, and furthermore, modes corresponding to differ¬ 
ent spin orientations are decoupled. Truncating the system ( |45| ) 
to 




/Qi 

0 

n '2 

0 ) 





0 


0 

D' 


ctriO 

c\{t) 


Q' 

0 

Qi 

0 


c\iO 

7 


lo 

a '2 

0 

Di j 




(46) 


with ^2 = ^2 cos T] yields, for the initial condition clj ( 0 ) = 1 
and cJi{Q) = 0 otherwise, the exact solution 


clj (f) = cos(Qif) €05(0.2^) ~ i sin(Qif) cos(f22f), (47a) 

cii(f) = 0, (47b) 

c}(f) = - sin(Qif) sin(n 2 ) “ icos(Qif) sin(Q 2 ) - (47c) 

c{(f) = 0. (47d) 

The probabilities |cl)j(f)p and |c}(f)p follow as 

|c!i(f)P = cos2(n2f)> (48a) 

|c}(f)P = sin^(n2f). (48b) 


Thus, the probabilities oscillate with the Rabi frequency 


, cos 77 

1' = 1 __1, 


(49) 


which is known as the Rabi frequency of the Kapitza-Dirac 
effect with linearly polarized light, modihed by the factor cos 77 . 

As a consequence of ( |49l l, a nonvanishing ellipticity de¬ 
creases the Rabi frequency. In the case of circular polarization 
it even becomes zero; that is, Kapitza-Dirac scattering is com¬ 
pletely suppressed. To test this analytical prediction we htted 
the probability to hnd the electron in the quantum state 1 / 7 ^ j as 
obtained from the numerical solution of the Dirac equation to 
the function cos^(waf) with the ht parameter The numeri¬ 
cal values for Ua are compared to Qj for different values of 77 
in Fig.|^ which shows a clear cos 77 dependency, as predicted 
by our analytical considerations. 

As demonstrated above, the dynamics of the electron in the 
antirotating setup can be modeled by a time-independent scalar 
potential, which may be understood by classical arguments. 
Consider a charged particle in the antirotating helds. The clas¬ 
sical motion of the particle is determined by the Lorentz force 


f = — + r X . (50) 

m 


The velocity r may be divided into two parts, r = + r^, 

that are parallel to the x axis and perpendicular to the x axis, 
respectively. The change of the latter is primarily determined 
by the electrical held 


= (51) 

m 

Integrating this equation of motion yields the perpendicular 
velocity 


t-± 


2qE 

77707 


sin(o7f) 


0 

cos{kx) 
^cos{kx + 77 )^ 


+ ri(0). 


(52) 


The motion along the x axis is solely determined by the mag¬ 
netic held, i. e.. 



FIG. 6: (Color online) Rabi oscillation frequencies for different de¬ 
grees of ellipticity 77 for the antirotating setup as obtained by a numer¬ 
ical ht procedure to data of numerical solutions of the Dirac equation 
compared to the predictions by the relativistic Pauli equation ( |17) . 
Parameters are as in Fig.[^ 
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= —/"i X t) 


4q^E^ sin^(wf) 


/ 0 ) 


( 0 

cos{kx) 

X 

- sin(kx + 77 ) 

\cos{kx + q)! 


^ sin(kx) y 


+ ^r,{0)xB^^{x,t)e,. (53) 
m 


Because the velocity rj^ is in phase with the magnetic field and 
both are proportional to sin(a;f), the acceleration parallel to 
the X axis does not average out over a laser cycle. Only the 
effect of the last term in Eq. ( |53| l cancels for sufficiently long 
interaction times. The value of the nonvanishing contribution 
depends crucially on the spatial phase relation between the 
electric field, which is ~ ( 0 ,cos(^jic),cos(A:x + and the 
magnetic field, which is ~ (0, - sin(^x + ri),sin{kx)y. Thus, 
the time-averaged acceleration parallel to the x axis 


2q^E^ 

r II Rj —--cos 77 sin( 2 ^x + 77 ) 


(54) 


and the corresponding potential 
2q^E^ 


V, 


pond ■ 


mco^ 


cos 77 cos‘^{kx + 77 / 2 ) 


(55) 


are proportional to cos 77 , which causes the cos 77 dependency 
of the Rabi frequency (|49|). 


5. Conclusions 

We studied the two-photon Kapitza-Dirac effect in the Bragg 
regime for setups with counterpropagating elliptically polar¬ 
ized electromagnetic waves whose electric-field components 
may have the same or an opposite sense of rotation. The time- 
dependent fully relativistic Dirac equation was solved numer¬ 
ically to simulate the electron’s dynamics. To allow for an 
analytical description of the scattering process the explicitly 
time-dependent Dirac-equation Hamiltonian was reduced to 
an effective time-independent eight-level Hamiltonian or four- 
level Hamiltonian by approximating the Dirac equation by the 
Pauli equation plus the leading relativistic corrections and em¬ 
ploying suitable ponderomotive potentials. In the parameter 
range of the Bragg regime, the fully relativistic Dirac equation 
and the effective time-independent Hamiltonian give qualita¬ 
tively the same results, which also agree quantitatively very 
well for not too long interaction times. 

The light waves’ ellipticity does not change the Bragg con¬ 
dition compared to the case of linear polarization. It induces, 
however, spin dynamics in the case of corotating counterprop¬ 
agating waves. The spin’s expectation value of the quantum- 
mechanical superposition of the scattered and nonscattered 
states precesses with a period which is larger than the Rabi 
period of the Kapitza-Dirac scattering. The projections of the 
electron’s quantum state onto the scattering and nonscatter¬ 
ing channels, however, show spin precession on the time scale 
of the Rabi period of Kapitza-Dirac scattering. The preces¬ 
sion of the electron’s spin is induced equally by the Zeeman 


interaction and the electron-spin’s coupling to the photonic 
spin density of the standing light wave. This effect may be 
observable by employing upcoming laser sources for circularly 
polarized light ll45l with wavelength and intensity parameters 
as utilized by Freimund et al. El, provided that the electron 
interacts long enough with the electromagnetic field to rotate 
the spin. Compared to nonresonant scattering Il29l[^ . where 
spin precession is also induced by the elliptically polarized 
light wave, the spin-precession frequency is substantially faster 
for the Kapitza-Dirac effect. Thus, shorter interaction times 
are required to observe spin-precession experimentally. 

Antirotating counterpropagating waves yield a standing wave 
with linear polarization, i. e., vanishing photonic spin density; 
thus, no spin effects can be observed for this kind of setup. 
However, the phase between the electric and the magnetic-field 
components of the standing wave depends on the ellipticity of 
the counterpropagating laser fields, and in this way the Rabi pe¬ 
riod of Kapitza-Dirac scattering depends on the ellipticity. This 
Rabi period diverges in the limit of circular polarization; that is, 
the Kapitza-Dirac effect cannot be observed in the antirotating 
setup for circular polarization. 


A. Nonresonant scattering 


Recently, it was shown m that the coupling of the spin angular 
momentum of light beams with elliptical polarization to the 
spin degree of freedom of free electrons can lead to spin preces¬ 
sion. The spin-precession frequency was derived by a tedious 
calculation based on time-dependent perturbation theory. The 
setup that was considered in Refs. Il29ll3()l is identical to the 
corotating field configuration, which wasO analyzed In Sec. |4.1[ 
except that the electron is initially at rest; that is, it does not 
fulfill the Bragg condition of the Kapitza-Dirac effect. In the 
following we will show that the spin-precession frequency of 
lOOll can be derived easily by employing the time-independent 
Hamiltonian in Eq. G3- 

Equation ( |2^ yields, for the even modes, the truncated sys¬ 
tem 


4(0 
4(0 
4(0 
\4(0 ) 


MQi 

0 

02 

n'j 

0 

0 

0 

4Di 


O 2 

0 

0 

Q 2 


0 

0 

02 

d;, 


0.2 

0 

0 


O 2 

0 

0 

O 2 


4Qi 

0 

1 0 

0 

0 '^ 

O 2 

0 

40 


' 42 ( 0 ' 

42(0 

4(0 

4(0 

4(0 

k 4(0 ) 


(A.l) 

As for the Kapitza-Dirac effect the eigenvalues and the eigen¬ 
vectors of the coefficient matrix in Eq. ( |A. 1[ ) can be calculated 
exactly, which gives the exact eigenvalues 


£ 1,2 = 2Di - ^jA^\ + 2{D.2±Q!^f , (A.2a) 

£ 3,4 = 4Qi, (A.2b) 

£ 5,6 = 2Di + + 2(02 ± ^ 3 )^ . (A.2c) 

However, we replace the exact eigenvectors by the approximate 
eigenvectors 
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V6) = 


/ 0 

0 

1 

n/2 

1 

0 


0 

0 

1 

^/2 

1 

■ n/2 

0 


\/2 

0 


^ 0 
1 

0 


0 0 


0 

0 

1 

■ n/2 

0 


0 

0 

1 

■ n/2 


The evolution of the probabilities |cj(f)p and |cg(f)p is with 
dAlli 


-5/ 

(A.3) 

Specifying cj(0) = 1 and c^(0) = 0 otherwise as the initial 
condition and approximating the eigenvectors of the coefficient 
matrix in Eq. ( |A.l| l by ( |A.3| ) yields the time-dependent solution 
todATll, 


|cj(f) p = cos^ ( 

= ■ (A-5b) 

Considering again the regime (|24ll, the approximation 


si-si 


2£l2£l'i 

Q, 


(A.6) 


holds. Thus, 


nJ —»■ 

II 

O 

(A.4a) 



02 ( 1 ) = 0 , 

(A.4b) 

1 t/ m2 2 ( q Sint; f \ 

ko(OI -cos Y)5;i2,„2c5 2 j ’ 

(A.7a) 


(A.4c) 

1 i/ m 2 . 2 / sin77 f \ 

ko(OI -sm (27r)5rr2m2c5 2 ) ’ 



(A.4d) 

(A.7b) 

c\{t) = 0, 

(A.4e) 

which represents exactly the Rabi oscillation of the electron 

0 

II 

—♦■rd 

0 

(A.4f) 

spin as predicted in Refs. 112911301. 
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